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Rate  constants  and  product  ion  branching  fractions  are  reported  for  the  reactions  of  CH3+,  C2Hs+,  1y-C3H7+, 
•y-C4H9+,  f-C4H9+,  and  f-CsHn*  with  O2  and  O3  at  300  K  in  a  variable-temperature  selected-ion  flow  tube 
(VT-SIFT).  The  reaction  rate  constant  for  CH3+  with  03  is  large  and  approximately  equal  to  the  thermal 
energy  capture  rate  constant  given  by  the  Su-Chesnavich  equation.  The  C2H5+,  j-C3H7+,  and  ^C4H9+  ions 
are  somewhat  less  reactive,  reacting  at  approximately  7—46%  of  the  thermal  capture  rate.  The  HCO+  and 
C2H30+  ions  are  the  major  products  in  these  reactions.  The  f-C4H9+  and  f-C5Hn+  ions  are  found  to  be  unreactive, 
with  rate  constants  <5  x  10-12  cm3  s”1,  which  is  the  present  detection  limit  of  our  apparatus  using  this  ozone 
source.  Ozone  is  a  singlet  in  its  ground  state,  and  ab  initio  calculations  at  the  B3LYP/6-31G(d)  level  of 
theory  indicate  that  reactant  complexes  can  be  formed,  decreasing  in  stability  with  the  size  of  alkyl  chains 
attached  to  the  cationic  carbon  atom.  The  decreasing  reactivity  of  the  alkyl  ions  with  increasing  order  of  the 
carbocation  is  attributed  to  a  greatly  reduced  03  binding  energy.  The  ions  listed  above  do  not  undergo  two- 
body  reactions  with  02,  k  <  5  x  10“13  cm3  s-1,  despite  the  availability  of  reaction  channels  with  exothermicities 
of  several  hundred  kilojoules  per  mole.  Ab  initio  calculations  at  the  B3LYP/6-31G(d)  level  of  theory  indicate 
that  the  02  reaction  systems  form  weak  complexes  with  large  C— O  bond  distances  (repulsive  at  smaller 
distances)  on  the  lowest  energy  triplet  potential  energy  surface.  Access  to  the  singlet  surface  is  required  for 
bond  formation;  however,  this  surface  is  not  accessible  at  thermal  energies. 


I.  Introduction 

Carbonium  ions  are  found  in  a  wide  variety  of  environments, 
including  interstellar  clouds,  chemical  vapor  deposition  plasmas, 
and  combustion  environments.  Studies  of  the  oxidation  of  these 
ions  have  historically  been  focused  on  CrtHm+  ion  reactions  with 
O  atoms.  These  reactions  have  been  studied  to  elucidate  the 
synthetic  chemistry  occurring  in  interstellar  clouds.  Fehsenfeld1 
examined  the  reaction  of  CH3+,  Bohme,  Mackay,  and  Schiff2 
examined  the  reaction  of  CHs+,  and  Viggiano  et  al.3’4  examined 
reactions  of  CH+,  CH5+,  C2+,  C2H+,  and  C2H2+  and  later 
isotopomers  of  CH3+  (CD„H3-n'f).  More  recently,  Le  Page  et 
al.5*6  reported  the  investigation  of  aromatic  hydrocarbon  ions 
Ci0H6+,  C10H7+  and  CioH8+  and  Ci6Hio+,  Ci6H9+,  and  Ci6Hn+ 
with  H,  N,  and  O  atoms.  Of  these  ions,  CioHg+,  Ci6Hio+,  and 
Ci6H9+  were  found  to  react  with  atomic  oxygen.  The  Ci6H9+ 
cation  showed  striking  reactivity  toward  atoms,  compared  to 
the  relative  reactivity  of  phenylium  and  naphthylium  ions,  which 
was  attributed  to  the  triplet  nature  of  the  Ci6H9+  ion  in  its  ground 
state.  Results  of  a  study  of  13  CmHrt+  ions  ( m  <  6)  with  O,  O2, 
and  NO  carried  out  by  Scott  et  al.7  showed  that  most  of  the 
hydrocarbon  ions  studied  exhibit  relatively  rapid  reactions  with 
O  atoms  that  proceed  at  substantial  fractions  of  the  collision 
rate.  One  notable  exception  found  by  Scott  et  al.  was  the  C2H5+ 
ion,  one  of  the  subject  ions  in  the  present  study,  which  showed 
no  reactivity  with  O  atoms. 

*  Corresponding  author. 

t  Permanent  address:  Department  of  Chemistry,  Montana  State  Univer¬ 
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The  formation  and  destruction  mechanisms  of  hydrocarbon 
ions  in  combustion  environments  is  an  area  of  recent  interest. 
In  particular,  plasma  sources  are  being  explored  as  ignition  and 
combustion  aides  in  supersonic  combustion,  where  fuel  activa¬ 
tion  and  initiation  on  fast-flow  time  scales  pose  a  daunting 
challenge.8  However,  the  information  gained  from  these  studies 
is  broadly  applicable  to  a  wide  range  of  problems  such  as  spark 
inhibition,  improved  engine  performance,  service  life,  explosion 
limits  in  blended  fuels,  hydrocarbon  molecular  growth,  and 
ignition.  The  oxidation  of  hydrocarbon  ions  in  these  environ¬ 
ments  is  of  fundamental  interest.  For  example,  C„Hm+,  C„HmO+, 
and  CrtHm02+  positive  ions  have  been  observed  in  the  flames 
of  a  variety  of  fuels.9  Although  these  ions  are  commonly 
observed,  their  formation  mechanisms  are  not  fully  understood. 

One  mechanism  involves  oxidation  of  hydrocarbon  ions  in  a 
single  step  via  reactions  with  oxygen  species  such  as  O,  02,  O2 
(a 1  A),  and  03.  As  mentioned  above,  some  work  has  been  carried 
out  with  regard  to  O  atom  reactions,  but  the  other  oxidation 
processes  remain  relatively  unexplored.  The  reactions  of  alkyl 
ions  with  numerous  compounds  have  been  reported  previously.10 
However,  reactions  leading  to  the  oxidation  of  CnH2*+i+  ( n  = 
2—5)  species  have  not  been  reported.  Oxidation  of  other  organic 
ions  with  03  has  been  studied  in  the  gas  phase  by  Mendes  et 
al.11  using  a  pentaquadrupole  mass  spectrometer.  The  collision 
energy  of  the  ion— molecule  reactions  studied  in  those  experi¬ 
ments  was  ca.  1  eV.  Mendes  et  al.11  found  that  neutral  ozone 
transfers  an  oxygen  atom  to  several  positive  ioiis,  namely,  radical 
cations  of  pryridines,  alkyl  halides,  and  halogen  cations.  A 
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Figure  1.  Overview  of  the  selected-ion  flow  tube  (SIFT)  with  ozonator  and  absorption  cell. 


considerable  amount  of  O2  was  present  in  these  experiments. 
However,  when  the  experiments  were  repeated  with  pure  O2, 
no  O-atom-transfer  products  were  observed,  confirming  that  the 
oxidation  occurs  from  O3  alone.  Ab  initio  calculations  were  also 
performed  as  part  of  the  work  by  Mendes  et  al.11  to  rationalize 
the  formation  of  the  reaction  products  observed.  An  interesting 
conclusion  of  that  work  is  that  simple  nucleophilic  addition  of 
ozone  followed  by  O2  loss  is  the  most  likely  mechanism  for 
O-atom  transfer  from  ozone  to  radical  cations  of  pryridines,  alkyl 
halides,  and  halogen  cations.  This  new  mechanistic  pathway 
was  suggested  to  be  an  important  pathway  for  the  generation 
of  chemically  interesting  but  difficult-to-generate  ionized  oxides. 

In  this  paper,  the  chemistry  of  ionized  hydrocarbons  reacting 
with  oxygen  and  ozone  is  investigated  at  300  K  using  fast  ion 
flow  tube  kinetics  instrumentation.  The  measurements  are 
conducted  with  a  selected-ion  flow  tube  (SIFT)  instrument  at 
ca.  0.5  Torr  helium  buffer  pressure,  and  the  kinetics  are  observed 
over  ca.  3—4  ms  reaction  time.  Specifically,  rate  constants  and 
product  ion  branching  fractions  are  reported  for  the  reactions 
of  CH3+,  C2H5+,  j- C3H7+  j-C4H9+,  f-C4H9+  and  f-C5H„+  with 
O2  and  03  at  room  temperature.  The  HCO+  ion  is  a  major 
product  in  many  of  these  reactions,  so  the  reactivity  of  HCO+ 
with  O2  and  O3  is  also  investigated.  Density  functional  calcula¬ 
tions  using  the  B3LYP  functional12-14  and  6-31G(d)  basis  set15 
are  also  performed  on  the  electronic  ground-state  potential 
surfaces  of  complexes  of  O2  and  O3  with  select  hydrocarbon 
radical  cations  to  elucidate  the  reaction  mechanism  leading  to 
the  oxidation  of  these  ions.  The  new  experimental  findings  are 
discussed  in  the  context  of  these  theoretical  calculations  and 
previous  work  where  available. 

II.  Method 

Experimental  Details.  The  measurements  were  made  at  300 
K  using  the  SIFT  instrument  at  the  Air  Force  Research 
Laboratory.  This  apparatus  is  described  in  detail  elsewhere,16*17 
and  the  modifications  required  to  perform  experiments  with 
ozone  as  a  reactant  gas  have  been  discussed  recently.18  The  SIFT 
is  shown  schematically,  with  modifications,  in  Figure  1. 

In  brief,  carbonium  ions  are  prepared  in  a  remote,  differen¬ 
tially  pumped  chamber  by  electron  impact  ionization  of  alkyl 
bromides  and  alkyl  chlorides  entrained  within  a  supersonic  argon 
expansion.  After  a  skimmer,  all  ions  are  directed  into  a 
quadrupole  mass  filter,  where  the  ionic  species  of  interest  is 
mass-selected.  The  mass-selected  ions  are  injected  via  a  Venturi- 
type  inlet  into  a  fast  flow  of  helium  carrier  gas  in  a  meter-long 
stainless  steel  flow  tube.  Reactant  gas  is  introduced  into  the 
flow  tube  through  either  of  two  stainless  steel  reactant  gas  inlets 
and  allowed  to  react  over  a  known  distance  at  a  known  flow 
velocity.  The  SIFT  operates  at  ca.  0.5  Torr  helium  buffer 
pressure,  and  the  kinetics  are  observed  over  ca.  3—4  ms  reaction 


time.  A  second  quadrupole  mass  spectrometer  resolves  the 
reactant  and  the  product  ions,  which  are  then  detected  by  an 
electron  multiplier.  Extrapolations  of  product  branching  fractions 
to  zero  reactant  flow  yield  the  nascent  branching  fraction.  The 
decay  in  the  primary  reactant  ion  signal  as  a  function  of 
increasing  reactant  gas  flow  rate  yields  the  reaction  rate  constant. 
Concentrations  are  such  that  [buffer]  »  [reactant  neutral]  » 
[ions].  Under  these  conditions,  pseudo-first-order  kinetics  apply, 
and  the  rate  constant  is  given  by 


=J_, 

[B]r  n[A±] 


(1) 


where  k  is  the  rate  constant,  r  is  the  reaction  time,  [B]  is  reactant 
neutral  concentration,  and  [Aq^]  and  [A^]  are  the  primary 
reactant  ion  concentrations  in  the  absence  and  in  the  presence 
of  reactant  neutral,  respectively.  The  reaction  time,  r,  is  the 
reaction  distance  divided  by  the  buffer  velocity  multiplied  by  a 
correction  factor  determined  from  previous  time-of-flight  mea¬ 
surements  that  accounts  for  the  fact  that  both  the  ion  velocity 
and  the  ion  concentration  are  at  a  maximum  along  the  axis  of 
the  flow  tube.  A  typical  value  for  the  correction  factor  is  1.6. 
The  buffer  velocity  is  obtained  from  the  mass  flow  rate  of  the 
buffer,  the  flow  tube  cross  section,  temperature,  and  pressure 
in  the  normal  manner.19  The  absolute  uncertainties  of  the  rate 
constants  are  25%,  and  relative  uncertainties  are  15%. 

An  Orec  O3V-0  ozonator  was  interfaced  to  the  SIFT 
apparatus  to  produce  the  ozone  reactant  gas  used  in  this 
experiment.  A  commercial  supply  of  02  (Airco,  99.999%)  was 
used  to  manufacture  the  ozone.  The  03  was  produced  at  a 
pressure  of  3-6  psig  using  a  0.9-A  discharge  current.  The 
resulting  reactant  gas  is  approximately  5%  03  in  O2,  which  is 
the  same  as  that  reported  previously  by  Fahey  et  al.20  The 
fraction  of  O3  was  relatively  independent  of  the  reactant  gas 
flow  rate  over  a  100  SCCM  flow  range.  The  presence  of  02 
does  not  affect  the  measurements,  because  O2  was  found  to  be 
unreactive  with  the  positive  ions  studied  here,  including  the 
product  ions.  Flow  contamination  was  less  than  0.05%,  typically 
arising  from  nitrogen  and  carbon  dioxide  trace  gases  in  the 
ozonator  and  the  tubing.  Flushing  the  lines  and  the  ozonator 
regularly  to  prevent  the  buildup  of  N02,  N03,  and  C02 
minimized  contamination.  All  fittings  and  valves  used  in  the 
O3  delivery  system  are  stainless  steel,  and  the  tubing  was  either 
stainless  steel  or  Teflon. 

A  flow  controller  regulates  the  03/02  reactant  mixture  flow 
into  a  10.2-cm-long,  1.3-cm-diameter  Pyrex  absorption  cell  fitted 
with  quartz  windows.  The  absorption  cell  is  connected  to  the 
flow  tube  stainless  steel  reactant  inlets  by  approximately  40  cm 
of  0.25-in.-o.d.  Teflon  tubing.  Varying  the  length  of  the  Teflon 
tubing  from  40  to  300  cm  produced  no  difference  (<5%)  in 
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the  results,  suggesting  that  no  significant  03  decomposition 
occurs  in  the  tubing.  The  absolute  concentration  of  03  is 
measured  by  optical  absorption  at  248  and  254  nm  using  a 
Perkin-Elmer  Lambda  10  UV/vis  spectrometer.  The  concentra¬ 
tion  of  the  ozone  in  the  flow  tube  was  found  using 


where  A  is  the  log  absorbance  (base  10)  output  from  the 
spectrometer,  o\  (cm2  molecule"1)  is  the  absorption  cross  section 
for  O3  at  wavelength  X ,  /  is  the  length  of  the  absorption  cell  in 
cm,  Pft  is  the  flow  tube  pressure,  Pac  is  the  absorption  cell 
pressure,  Fac  is  the  total  flow  through  the  absorption  cell,  and 
Fft  is  the  total  flow  through  the  flow  tube.  The  248-  and  254- 
nm  cross  sections  used  were  1.08  x  10-17  and  1.137  x  10" 17 
cm2  molecule-1,  respectively.21*22  In  our  previous  publication,18 
results  regarding  several  negative-ion  reactions  with  ozone  were 
compared  to  literature  values  regarding  the  negative-ion  chem¬ 
istry  of  ozone,  which  prior  to  our  study  lacked  sufficient 
agreement  for  application  to  detailed  physical  models.  For  the 
reaction  of  SF6-  with  ozone,  there  is  now  good  agreement 
between  three  recent  measurements,18  indicating  reliable  delivery 
and  knowledge  of  the  O3  concentrations  in  the  present  experi¬ 
ment. 

CH3+  was  prepared  by  ionizing  methyl  bromide  (Aldrich, 
>99.5%),  C2H5+  was  prepared  by  ionizing  ethyl  bromide 
(Aldrich,  >99%),  C3H74"  was  prepared  by  ionizing  isopropyl 
bromide  (Aldrich,  >99%),  W-C3H7+  was  prepared  by  ionizing 
n-propyl  bromide  (Aldrich,  >99%),  t-C 4H9+  was  prepared  by 
ionizing  terf-butyl  chloride  (Aldrich,  >99%),  j-C4H9+  was 
prepared  by  ionizing  n-butyl  chloride  (Aldrich,  >99%),  and 
C5Hn+  was  prepared  by  ionizing  2-chloro-2-methylbutane. 
HCO+  was  formed  by  ionizing  CO  (Airco,  99.9%)  in  the 
supersonic  source  and  injecting  CO+  into  the  flow  tube,  which 
is  reacted  with  10  SCCM  of  H2  added  20  cm  upstream  of  the 
reactant  inlet.  All  source  gases  were  entrained  within  a 
supersonic  argon  (Airco,  99.995%)  expansion. 

Calculations.  Geometry  optimizations  were  carried  out  using 
Gaussian’s  implementation  of  the  B3LYP  hybrid  density 
functional12-14  using  the  standard  6-31G(d)  basis  set  of  Pople 
et  al.15  Open-shell  fragments  were  treated  using  the  spin- 
unrestricted  density  funcional  formalism.  Analytical  frequency 
calculations  were  used  to  characterize  the  nature  of  these 
stationary  points.  All  of  the  reaction  pathways  were  verified 
by  intrinsic  reaction  coordinate  (IRC)  calculations23-27  both 
forward  and  backward  from  the  encountered  transition  states. 
Standard  Mulliken  population  analysis28-30  was  employed  to 
analyze  electron  density  and  spin  density  distributions  from  the 
B3LYP/6-31G(d)  densities. 

III.  Results 

Reaction  rate  constants  and  product  branching  fractions  for 
alkyl  ions,  CJEWi*  ( n  =  1-5),  and  HCO+  reacting  with  03 
measured  at  300  K  with  the  SIFT  are  shown  in  Table  1 .  None 
of  the  subject  cations  showed  any  measurable  bimolecular 
reactivity  with  02,  k  <  5  x  10-13  cm3  s-1.  In  Table  1,  the 
standard  reaction  enthalpies  at  298  K  (in  kJ  mol-1)  have  been 
calculated  and  are  listed.  The  heats  of  formation  are  taken  from 
the  NIST  compilation.31,32  Note  that  the  terr-butyl  ion  heat  of 
formation  has  been  recently  revised,  and  a  value  of  A//f°298(r- 
C4H94")  =  711  kJ/mol  from  recent  experimental33-35  and 
theoretical  studies36  is  used,  which  is  20  kJ/mol  higher  than 
previous  experimental  measurements.  Figure  2  shows  corrected 


03  Reactant  Concentration  (xIO11  cm"3) 


Figure  2.  Corrected  (see  text)  primary  ion  decay  plots  for  CH3+  (■), 
HCO+  (A),  C2H5+  (•),  S-C3ITF  (♦),  and  r-C4H9+  (*)  reacting  with 
O3.  The  solid  lines  are  nonlinear  least-squares  fits  to  the  data,  performed 
to  determine  the  rate  constants  according  to  eq  1.  Experiments  are 
performed  with  a  reactant  gas  composition  of  5%  03  in  02,  which  is 
possible  since  02  is  unreactive. 


03  Reactant  Concentration  (xIO12  cm"3) 


Figure  3.  Normalized  counts  (■)  plotted  vs  O3  concentration. 
The  dashed  line  is  the  result  of  a  single-exponential  fit  of  the  data  to 
determine  the  rate  constant  according  to  eq  1.  The  result  of  the  single¬ 
exponential  fit  is  labeled  Xrtotai-  The  solid  line  is  the  result  of  a  double¬ 
exponential  fit  of  the  data  which  yields  a  fast  component  ( ksec )  with  a 
rate  constant  of  (6  ±  4)  x  10-n  cm3  s-1  and  a  slow  component  (Jfctert) 
with  a  rate  constant  <5  x  10" 12  cm3  s"1. 

(see  below)  primary  ion  decay  plots  for  the  reactions  of  CH3+, 
HCO+,  C2Hs+,  ^-C3H7+,  and  r-C4H9+  with  03.  The  solid  lines 
are  linear  least-squares  fits  to  the  data,  performed  to  determine 
the  rate  constants  according  to  eq  1.  Sample  kinetics  data  for 
C4H9+,  prepared  by  ionizing  n-butyl  bromide,  are  shown  in 
Figure  3.  Electron  impact  ionization  of  n-butyl  bromide  is 
expected  to  produce  both  s-C4Il9+  and  M^Hg"1".37  The  isomeric 
composition  was  determined  by  exploiting  the  order  of  mag¬ 
nitude  difference  in  the  hydride-transfer  reaction  rate  constant 
with  2-methylpentane  (i-CsH^)  for  the  two  isomers,  as  shown 
in  Figure  4.  The  specific  results  for  each  cation  studied  are 
discussed  below. 

HCO4*  Cation.  The  HCO+  cation  is  a  major  product  of  the 
primary  and  secondary  carbocation  reactions  with  03  listed  in 
Table  1.  Hence,  the  reactivity  of  HCO+  with  O3  was  included 
in  this  study.  In  these  studies,  CO+  was  injected  into  the  flow 
tube  and  reacted  with  10  SCCM  of  H2  added  at  inlet  1.  This 
provided  a  source  of  HCO+,  although  other  ions  including  H30+ 
and  02+  were  present  in  trace  amounts.  The  SIFT  spectrum  of 
CO+  in  the  absence  of  added  H2  indicated  the  presence  of  H20+ 
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TABLE  1:  Reaction  Rate  Constants  for  Reactions  of  Ozone  at  300  K,  Measured  with  the  Selected-Ion  Flow  Tube  (SIFT) 

k,[kc] 

branching 

reaction 

products 

(x  10  9  cm3  s-1) 

fractions 

AH  (kJ/mol) 

hco*  +  o3— 

H<v  +  CO 

1.2,  [U] 

>0.98 

-31.5 

ch3++o3— 

HCO+  +  (H2  +  02) 

1.7,116] 

0.66 

-410 

h2co+  +  ho2 

0.16 

-293 

H202+  +  HCO 

0.10 

-311 

o2+  +  ch3o 

0.07 

-54 

h3o+  +  co2 

<0.01 

-1038 

C2Hs++03  — 

HCO+  +  (CH4  +  02) 

0.59,  [1.3] 

0.70 

-294 

C2H30+  +  (H2+02) 

0.30 

-392 

j-C3H7+  +  03— 

C2H30+  +  (CH„  +  02) 

0.22,  [1.2] 

0.64 

-363 

HCO+  +  (C2H6  +  02) 

0.16 

-200 

CH30+  +  C2H4  +  02 

}o.l4 

-186 

CH30+  +  CHjCOOH 

-670 

C3H30+  +  2H20 

0.05 

-674 

c3h6o+  +  ho2 

0.01 

-167 

s-C4H9+  +  03— 

HCO+  +  (C3H8  +  02) 

0.08,  [1.1] 

major  (~0.40) 

-188 

c2h3o+  +  (C2H6  +  02) 

major  (~0.35) 

-339 

c2h5o+  +  c2h4  +  o2 

-273 

C2H50+  +  ch3cooh 

jminor  (~0.14) 

-758 

c2h2o++c2h6  +  ho2 
ch3o+,c3h3o+,  c3h5o+ 

minor  (~0.09) 
trace  (<0.02) 

-111 

/-C4fV+03- 

C2H30+  +  (C2H«  +  02) 

<0.005,  [1.1] 

-284 

HCO+  +  (C3H8  +  02) 

-133 

f-C5H,|++  03  — 

C2H5C0+  +  (C2H6+02) 

<0.005,  [1.1] 

-296 

c2h3o+  +  (C3H8  +  02) 

-255 

HCO+  +  (C4H10  +  02) 

-112 

All  of  the  cations  listed  showed  no  reactivity,  k<  5  x  10“13  cm3  s_1,  with  02.  The  calculated  collision  rate  constant,  kc,  and  the  measured  rate 
constant,  k ,  are  listed  in  italics.  The  reaction  products  and  branching  fractions  are  listed  under  the  rate  constants.  Energetics  are  taken  from  the 
NIST  Chemistry  WebBook.31*32  Neutral  products  are  not  resolved  in  these  experiments.  Brackets  indicate  that  more  than  one  product  channel  could 
contribute  to  the  branching  fraction.  Products  contained  in  parentheses  are  estimated  on  the  basis  of  the  mechanism  discussed  in  the  text  involving 
02  loss.  Note  that  the  reactions  are  considerably  more  exothermic  if  the  reaction  products  in  parentheses  are  not  dissociated. 


Figure  4.  Normalized  C4IV  counts  (■)  plotted  vs  2-methylpentane 
(/-C5H12)  concentration.  The  solid  line  is  the  result  of  a  double¬ 
exponential  fit  of  the  data  which  yields  a  fast  component  (£**)  with  a 
rate  constant  of  (4.9  ±  1.0)  x  10~10  cm3  s-1  and  a  slow  component 
(fcten)  with  a  rate  constant  (4.2  ±  0.3)  x  10"11  cm3  s"1. 

and  02+,  which  indicates  that  the  source  chemistry  is  responsible 
for  the  background  levels  of  H30+  and  C>2+.  HCO+  is  known 
to  be  unreactive  with  O2,38  and  that  finding  was  confirmed  in 
this  study.  However,  HCO+  is  very  reactive  with  O3,  which  is 
expected  since  the  isoelectronic  hydrogen  cyanide  molecule  is 
known  to  react  with  ozone  to  produce  HCNO  as  a  major 
product.39  The  major  product  of  HCO+  reacting  with  O3  is 
H03+.  The  branching  fraction  reported  in  Table  1  is  >0.98  and 


is  intended  to  reflect  the  uncertainty  of  the  source  conditions. 
The  measured  rate  constant  of  1.2  x  10-9  cm3  s-1  is  nearly 
equal  to  the  thermal  capture  rate  constant  of  1.3  x  10“9  cm3 
s'1. 

Methyl  Cation.  The  methyl  cation  was  prepared  by  ionization 
of  methyl  bromide  and  was  selectively  injected  into  the  flow 
tube.  However,  despite  very  low  injection  energies,  approxi¬ 
mately  12%  of  the  CH3+  collisionally  dissociated  to  form  CH+ 
and  CH2+.  While  these  breakup  products  do  not  interfere  with 


the  reaction  rate  constant  determination,  they  do  react  with  02 
to  form  some  of  the  same  products  observed  in  the  CH3+  +  03 
reaction.40  Namely,  CH+  and  CH2+  react  with  02  with  rate 
constants  of  9.7  x  10“ 10  and  9.1  x  10~"10  cm3  s”1,  respectively, 
and  produce  the  following  products:40 

CH+  +  02  —  HCO+  +  O 

(major) 

(3a) 

—  0+H-HC0 

(minor) 

(3b) 

—  CO+  +  OH 

(minor) 

(3c) 

and 

CH2+  +  02  —  HCO+  +  OH 

(major) 

(4a) 

-*h2co+  +  o 

(minor) 

(4b) 

The  methyl  cation,  on  the  other  hand,  does  not  undergo  two- 
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body  reactions  with  O2  and  only  exhibits  a  slow  ternary 
association  reaction,  &300K  “  8.6  x  10“30  cm6  s"1. 40,41 

Therefore,  the  impurity  CH+  and  CH2+  ions  were  eliminated 
before  entering  the  reaction  zone  by  addition  of  a  small  amount 
(5  SCCM)  of  O2  approximately  20  cm  upstream  of  the  first 
reactant  inlet.  The  0+  and  CO+  ions  resulting  from  reactions 
3b  and  3c  are  converted  to  02+  via  charge  transfer  with  02, 
and  the  H2CO+  from  reaction  4b  is  converted  to  HCO+  (major) 
and  H2O24"  (minor)  via  reaction  with  O2.  After  the  addition  of 
O2  at  inlet  1,  the  resulting  ion  composition  in  the  flow  tube  at 
the  reactant  inlet  was  observed  to  be  approximately  88%  CH3+, 
8.5%  HCO+,  2.5%  02+,  and  -1%  H202+.  The  02+  and  H202+ 
ions  were  accounted  for  by  subtracting  their  baseline  intensities 
observed  at  zero  02/C>3  reactant  flow  from  the  appropriate 
product  ion  intensities.  The  HCO+  ion  counts  observed  at  zero 
02/03  reactant  flow  were  corrected  for  by  subtracting  the 
appropriate  number  of  counts  at  each  ozone  concentration,  on 
the  basis  of  the  reactivity  of  HCO+  with  O3  listed  in  Table  1. 

The  methyl  cation  reacts  at  the  collision  rate  with  O3, 
producing  HCO+  as  a  major  product  and  H2CO+,  H202+,  02+, 
and  H30+  as  minor  products.  The  product  ion  counts  observed 
for  H2CO+,  H202+,  and  02+  were  significantly  above  the 
background  counts  discussed  for  these  ions.  Therefore,  these 
ions  are  confirmed  products.  The  H30+  product,  on  the  other 
hand,  results  either  from  direct  reaction  of  CH3+  with  O3,  which 
is  highly  exothermic  but  involves  significant  rearrangement,  or 
from  secondary  reactions  of  HCO+  with  the  very  low  level  of 
H20  present  in  the  flow  tube.  Secondary  chemistry  can  be 
accounted  for  by  extrapolating  the  baseline-corrected  product 
branching  fraction  back  to  zero  reactant  concentration.18  How¬ 
ever,  this  approach  did  not  fully  account  for  the  amount  of  H30+ 
observed.  Therefore,  H30+  is  listed  as  a  trace  product  in  the 
reaction  of  CH3+  +  O3. 

Ethyl  Cation.  The  ethyl  cation,  prepared  by  ionizing  ethyl 
bromide,  reacts  with  O3  to  form  two  primary  products,  HCO+ 
and  C2H30+.  Because  C2Hs+  and  HCO+  are  mass-coincident, 
two  separate  experiments  were  performed  to  determine  the 
reaction  rate  constant.  One  experiment  involved  the  addition 
of  CH3CI  to  the  flow  tube  via  inlet  1.  Methyl  chloride  reacts 
with  HCO+  by  proton  transfer  to  form  CEUC1+  and  was  not 
intended  to  react  with  C2H5+.  However,  the  addition  of  methyl 
chloride  did  cause  a  reduction  in  the  C2H5+  signal  and  produced 
several  unidentified  high-mass  ions  containing  chlorine,  so  the 
amount  of  methyl  chloride  added  was  such  that  only  about  20% 
of  the  original  C2H5+  signal  remained.  In  a  second  experiment, 
13CCH5+  ion  was  injected  into  the  flow  tube.  Statistically  equal 
amounts  of  12CHO+  and  I3CHO+  are  assumed  to  be  formed, 
and  the  intensity  of  12CHO+  is  used  as  a  measure  of  the  13CHO+ 
contribution  to  the  mlz  =  30  (sum  of  13CCH5+  +  13CHO+). 
The  reaction  rate  constant  was  determined  by  subtracting  the 
intensity  monitored  at  mlz  29  (12CHO+)  from  that  at  mfz  —  30 
(sum  of  13CCH5+  4-  13CHO"1')  and  plotting  the  decay  of  the 
resulting  ion  signal  as  a  function  of  O3  flow.  The  two  different 
experiments  produced  the  following  results.  Using  the  methyl 
chloride  scavenger,  three  determinations  yielded  rate  constants 
of  6.3  x  10”10,  6.5  x  10”10,  and  5.7  x  10“10  cm3  s”1.  The 
^CCHs4  experiment  in  a  single  determination  yielded  5.2  x 
10“ 10  cm3  s”1.  The  reaction  rate  constant  in  Table  1  is  the 
average  of  the  four  measurements. 

The  reaction  product  branching  fraction  was  estimated  by 
reacting  the  12C  form  of  C2H5+  with  high  O3  concentrations 
where  no  intensity  at  mlz  =  29  remained,  i.e„  C2Hs+  totally 
reacted  and  all  of  the  HCO+  was  converted  to  H03+.  Ap¬ 
proximately  5%  of  the  C2H5+  dissociated  to  C2H3+  upon 


injection  into  the  flow  tube.  The  reaction  products  originating 
from  C2H3+  were  not  identified  or  corrected  for  and  constitute 
no  more  than  a  5  percentage  point  uncertainty,  i.e.,  the  original 
fraction  of  the  C2H3+  reactant  ion  intensity,  assuming  it  yields 
only  a  single  product  in  the  reported  values.  Since  there  are  no 
minor  channels,  the  presence  of  C2H3+  has  a  small  effect  on 
the  overall  reported  branching  fractions. 

sec-Propyl  Cation.  The  propyl  radical  cation  was  prepared 
by  ionizing  isopropyl  bromide  in  the  supersonic  source.  Ap¬ 
proximately  5%  of  the  C3H7+  dissociated  to  C3H5+  upon 
injection  into  the  flow  tube.  Because  C2H30+  is  a  major  product 
of  the  jcc-propyl  cation  reaction  with  O3  and  is  mass-coincident 
with  C3H7+  the  reaction  rate  constant  was  determined  by 
injecting  the  13CC2H7+  (mlz  =  44)  form  of  the  reactant  ion. 
Given  that  the  expected  C2H30+  product  distribution  is  33% 
C2H30+  (mlz  =  43)  and  67%  13CCH30+  (mlz  =  44),  then  the 
/(13CC2H7+)  =  l(mlz  —  44)  —  2 1(mJz  =  43).  The  reaction  rate 
constant  was  determined  using  this  corrected  C3H7+  intensity. 
The  product  branching  fractions  were  determined  from  the  same 
13C  experiments  as  follows:  7(C2H30+)  =  3/(12C2HaO+), 
/(HCO+)  =  7(12CHO+)  +  7(13CHO+)>  and  7(CH30+)  =  7(12- 
CH30+)  +  7(i3CH30+),  with  the  remaining  C3  products  exhibit¬ 
ing  only  a  single  peak.  Three  determinations  were  averaged  and 
yielded  16%  CHO+,  14%  CH3O+,  64%  C2H30+,  5%  C3H3O+ 
and  1%  C3H60+.  The  reaction  products  originating  from  C3H5+ 
have  not  been  positively  identified  but  can  be  deduced  to  be 
either  CHO+  and/or  C2H30+,  on  the  basis  of  the  observed 
product  distributions  for  the  C3H7+  and  j-C4H9+  (see  below) 
reactant  ions.  Note  that  both  C3H7+  and  $-C4H9+  show  a  small 
amount  of  dissociation  to  C3H5+  upon  injection  into  the  flow 
tube  and  that  only  products  common  to  both  reaction  systems 
can  be  possible  reaction  products  of  C3H5+  and  O3.  The 
influence  of  the  C3H5+  on  the  reported  branching  fractions  has 
not  been  corrected  for  but  constitutes  no  more  than  a  5 
percentage  point  uncertainty  (original  fraction  of  the  reactant 
ion  intensity  yielding  a  single  product)  in  the  reported  values 
of  the  HCO+  and  C2H30'1'  product  channels. 

Previous  theoretical  studies  at  the  Hartree-Fock  and  MP2 
levels  of  theory  on  the  character  of  the  potential  energy  surface 
of  the  propyl  cation42  conclude  that  there  exist  two  global 
minima,  that  for  the  jec-propyl  cation  and  that  for  comer- 
protonated  cyclopropane.  These  calculations  showed  that  a 
minimum  does  not  exist  that  corresponds  to  the  n-propyl  cation, 
which  appears  to  be  only  a  transition  structure  in  the  intercon¬ 
version  of  the  sec-propyl  cation  and  the  comer-protonated 
cyclopropane  cation.  According  to  these  calculations,  the 
n-propyl  cation  is  80.5  kJ/mol  higher  and  the  comer-protonated 
cyclopropane  cation  is  30  kJ/mol  higher  in  energy  than  the  sec- 
propyl  cation  structure.  Our  B3LYP/6-31G(d)  calculations  are 
consistent  with  these  findings.  Furthermore,  experimental 
investigations  of  the  interconversion  of  the  propyl  cation 
structures  demonstrate  rapid  rearrangement  to  the  sec-propyl 
cation  form.  A  detailed  study  by  Ausloos  and  co-workers43 
showed  that  n-propyl  cations  isomerize  intramolecularly  to  either 
isopropyl  ions  or  protonated  cyclopropane  ions  within  10” 10  s. 
Rearrangement  to  the  isopropyl  ion  is  favored  and  increases  in 
importance  with  increasing  internal  energy  content  of  the  ion. 
Mass  spectrometric  investigations  suggest  that  the  isomerization 
of  the  comer-protonated  cyclopropane  structure  to  the  isopropyl 
structure  requires  10-7 — 10“ 5  s.44,45  All  of  these  isomerization 
times  are  much  shorter  than  the  approximately  1  ms  it  takes 
the  ions  to  leave  the  source  region  and  travel  to  the  reactant 
inlet.  Therefore,  the  reactivity  shown  in  Table  1  for  the  propyl 
cation  has  been  assigned  to  the  j-C3H7+  isomer. 
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To  confirm  this  hypothesis,  the  propyl  cation  was  also 
prepared  by  ionizing  n-propyl  bromide  in  the  supersonic  source. 
The  kinetics  measurements  were  performed  using  the  13CC2H7+ 
ion  and  corrected  for  13CCH30+  product  as  described  in  the 
experimental  details  using  isopropyl  bromide  in  the  source.  The 
product  branching  distribution  was  determined  from  the  mass 
spectrum  recorded  with  100  SCCM  of  O3,  where  no  C3H7+ 
remains.  As  above,  no  correction  for  the  products  originating 
from  C3H5+  (~5%  of  the  original  reactant  ion  intensity)  was 
made.  The  reaction  rate  constants  and  the  product  distributions 
for  the  C3H7+  ions  generated  from  n-propyl  bromide  are  found 
to  be  nearly  identical  to  those  produced  from  isopropyl  bromide. 
The  minor  differences  observed  in  the  branching  ratios  can  be 
attributed  to  statistical  deviations  in  the  reproducibility  of  the 
measurements  and  small  differences  in  the  source  and  injection 
conditions.  This  result  is  consistent  with  the  theoretical  and 
experimental  studies  discussed  above  regarding  the  stability  and 
isomerization  processes  of  the  C3H7+  ions,  which  indicate  that 
only  the  s-C3H7+  isomeric  form  of  the  ion  is  present  in  our 
experiments. 

tert- Butyl  Cation.  In  pulsed  ion  cyclotron  resonance  (ICR) 
experiments,  Shold  and  Ausloos37  observed  that  CJrLg*  cations 
formed  by  electron  impact  ionization  of  isobutane,  neopentane, 
2,2-dimethylbutane,  isobutyl  halides,  and  tert-butyl  halides  all 
have  the  tertiary  structure.  Fragmentation  of  n-alkanes,  2-me- 
thylbutane,  3-methylpentane,  n-butyl  halides,  and  sec-butyl 
halides  produce  both  s-CtHg*  and  r-C^*,  with  the  s-C4H9+ 
surviving  without  rearrangement  for  at  least  0.1  s.  However,  in 
the  case  of  the  halides,  a  collision-induced  isomerization  of  the 
j-C4H9+  to  the  r-C4H9+  was  found  to  occur.  The  ICR  experi¬ 
ments  of  Shold  and  Ausloos  were  conducted  at  pressures  of  ca. 
10“6  Torr  and  observation  times  ranging  from  10“3  to  0.5  s. 
The  temperature  in  the  analyzer  cell  was  kept  at  320  K,  and 
the  electron  energy  was  varied  between  10  and  25  eV.37 

On  the  basis  of  the  Shold  Ausloos  results,37  the  terf-butyl 
cation  was  prepared  by  electron  impact  ionization  of  tert- butyl 
chloride  in  the  supersonic  argon  expansion.  The  reported 
reaction  rate  constant  was  measured  by  increasing  the  flow  tube 
pressure  to  0.64  Torr  (throttling  the  roots  pump)  and  provided 
a  measured  value  of  <3.4  x  10“12  cm3  s"1,  which  is  at  the 
detection  limit  of  our  experiments  and  therefore  represents  an 
upper  limit.  No  reactivity  of  f-C4H9+  with  O2  was  observed 
either,  i.e.,  k  <  3  x  10~13  cm3  s_I. 

sec-Butyl  Cation.  The  sec-butyl  cation  was  prepared  by 
ionizing  n-butyl  bromide  in  the  supersonic  source.  Injection  of 
C4H9+  into  the  flow  tube  yielded  approximately  93%  C4H9+, 
5%  C3Hs+,  and  2%  C2H5+.  The  electron  impact  ionization  of 
n-butyl  bromide  produces  both  s-C4H9+  and  r-C4H9+,  and 
collisions  with  halide  molecules  in  the  source  region  can  convert 
^-C4H9+  to  /-C4H9+.37  The  kinetics  plots  with  the  G|H9+  cation 
exhibited  curvature  indicating  the  presence  of  more  than  one 
form  of  C4H9+.  Using  a  large  number  of  flow  points  and  fitting 
the  data  as  a  double  exponential  with  a  nonlinear  least-squares 
analysis  program  provides  an  estimate  of  the  relative  amounts 
of  the  two  types  of  C4H9+  and  a  measure  of  the  rate  constant 
for  the  fast-reacting  species.  Figure  3  shows  the  normalized 
counts  of  the  C4H9+  reactant  ion  plotted  as  a  function  of  the 
ozone  reactant  concentration.  The  result  of  a  single-exponential 
fit  to  the  data  yielded  the  dashed  line  associated  with  the  rate 
constant  designated  /:tota].  The  single-exponential  fit  reproduces 
the  data  very  poorly.  However,  a  double-exponential  fit  provides 
a  much  better  representation  of  the  data.  Free  fitting  several 
sets  of  data  for  the  two  rate  constants  yields  a  fast  component 
with  a  rate  constant  of  8  x  10“ 11  cm3  s'1  and  a  slow  component 


with  a  rate  constant  <4  x  1(T12  cm3  s'-1.  No  reactivity  of  either 
form  of  C4H9+  with  O2  was  observed,  i.e.,  k  <  3  x  10“13  cm3 
s”1,  which  is  the  detection  limit  of  the  experiment.  It  is  assumed 
that  the  slower  reacting  component  is  the  ferf-butyl  cation  (ca. 
67%)  and  the  faster  reacting  component  is  sec-butyl  cation  (ca. 
33%). 

Additional  insight  into  the  relative  amounts  of  the  different 
forms  of  C4H9+  was  explored  using  2-methylbutane  as  the 
neutral  reactant  in  hydride-transfer  reactions,  as  done  by  Shold 
and  Ausloos37  and  Meot-Ner  and  Field.46  The  tertiary  and 
secondary  forms  of  C4H9+  have  markedly  different  rates  for 
hydride  transfer  from  2-methylbutane  (Z-C5H12).  On  the  basis 
of  previous  results,10,46  the  hydride-transfer  rate  constants  for 
the  reactions 

f-C4H9+  +  (CH3)2CHCH2CH3  —  i-C4rf10  +  f-C5Hu+, 

k29tK  =  2  x  1(TU  cm3  s-  (5a) 

and 

j-C4H9+  +  (CH3)2CHCH2CH3  —  n-C4H10  +  r-C5H„+, 

k29S  k  =  ^  x  10“ 10  cm3  s'1  (5b) 

differ  by  over  an  order  of  magnitude.  Figure  4  shows  a  plot  of 
the  total  C4H9+  reaction  ion  concentration  as  a  function  of 
2-methylbutane  reactant  concentration.  A  biexponential  fit  of 
the  data  plotted  in  Figure  4  yields  two  rate  constants  of  (4.9  dr 
1.9)  x  1(T10  (35%)  and  (4.2  dr  1.4)  x  10“u  cm3  s^1  (65%). 
These  hydride-transfer  rate  constants  are  in  good  agreement  with 
those  reported  previously,  and  this  experiment  confirms  that 
approximately  35%  of  the  ions  formed  under  our  conditions 
are  S-C4H9+ 

The  product  branching  measurements  are  considerably  more 
difficult  than  those  of  the  reaction  rate  constant  because  most 
of  the  C4H9”**  is  the  nonreactive  tertiary  form.  Hence,  the 
contributions  of  contaminant  ions  produced  on  breakup  during 
injection,  C2H5+  and  C3Hs+,  are  considerably  larger  compared 
to  the  products  of  the  s-C4H9+  primary  ion.  The  major  products 
observed  are  HCO+  and  C2H30+,  and  C2H5O4*  and  C2H20+ 
are  minor  products.  As  shown  in  Table  1,  the  reaction  of  C2H54" 
with  03  produces  HCO+  and  C2H30+,  and  C2Hs+  is  unreactive 
with  02.  Since  C2H20+  and  C2H50+  were  not  observed  in  the 
experiments  discussed  above  that  contain  the  C3Hs+  impurity 
ion,  these  products  most  likely  originate  from  the  reaction  of 
s-C 4H94"  with  03.  Given  this  level  of  uncertainty,  the  only 
conclusions  that  can  be  drawn  are  that  HCO+  and  C2H30+  are 
major  products  and  that  C2H20+  and  C2H50+  are  minor 
products  produced  in  the  reaction  of  s-C4H9+  with  O3. 

The  production  of  C3H50+  in  the  reaction  of  s-C4H9+  with 
03  is  exothermic  and  is  a  conceivable  product  channel.  The 
reaction  for  such  a  process  is  given  below: 

j-C4H9+  +  O3—  c3h5o+  +  ch4  +  o2, 

AH  =  —393  kJ/mol  (6) 

Since  C3H50+  is  mass-coincident  with  C4H9+,  the  13CC3H9+ 
(mlz  ~  58)  form  of  the  reactant  ion  was  injected  into  the  flow 
tube.  Assuming  a  statistical  13C  product  distribution,  ap¬ 
proximately  25%  of  the  C3H50+  product  ion  produced  should 
be  C3H50+  (mlz  =  57).  Product  branching  fraction  measure¬ 
ments  taken  at  high  resolution  resulted  in  barely  detectable 
amounts  of  C3H30+  (mlz  =  55),  13C3H30+  (mlz  =  56),  and 
C3H50+  (mlz  —  57)  produced  in  the  reaction.  Trace  amounts 
of  ions  at  m/z  =  31  and  mlz  =  32  were  also  observed  and  have 
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been  assigned  to  CH30+  (m/z  =  31)  and  13CH30+  (mlz  =  32). 
Note  that  the  CH30+,  Cs^O4*,  and  C3HsO+  product  ions  could 
also  originate  from  reaction  of  the  C3Hs+  breakup  ion  (5%) 
reacting  with  03.  Therefore,  the  CH30+,  C3H30+,  and  C3H50+ 
product  ions  are  listed  as  trace  product  ions  in  Table  1 . 

C5H11  Cation.  The  CsHn+  radical  cation  was  prepared  by 
ionizing  2-chloro-2-methylbutane.  The  reported  reaction  rate 
constant  for  the  CsHn*  ion  formed  after  reaction  with  O2  is 
<3  x  10~13  cm3  s”1  and  with  03  is  <4  x  1(T12  cm3  s”1.  The 
tertiary  form  of  the  OsHa4-  cation  is  the  most  stable  form  of 
this  carbocation,  and  electron  impact  ionization  of  the  2-chloro- 
2-methylbutane  ion  precursor  is  expected  to  produce  this  ion 
exclusively.  The  r-CsHu*4  cation  is  the  most  stable  cation  studied 
in  this  series,47  and,  like  the  exceptionally  stable  cation, 

it  is  not  oxidized  by  either  O2  or  03,  even  though  several  very 
exothermic  products  are  available. 

IV.  Discussion 

Cacace  and  co-workers  first  experimentally  detected  the 
elusive  protonated  ozone  molecule,  H03+,  in  1994  by  reacting 
ozone  with  CHs+.48  The  present  study  shows  that  protonated 
ozone  is  produced  efficiently  by  the  reaction  of  HCO4*  with  03 
via  a  fast  proton-transfer  reaction,  as  listed  in  Table  1 .  The  high 
efficiency  is  expected  because  the  reaction  is  exothermic  and 
spin-allowed,  with  singlet  reactants  producing  singlet  products. 
There  are  several  other  examples  of  efficient  exothermic  proton- 
transfer  reactions  in  the  literature10  which  have  been  discussed 
in  detail. 

In  the  early  stages  of  combustion,  O2  is  in  large  abundance, 
and  even  limited  reactivity  of  hydrocarbon  ions  with  O2  could 
prove  to  be  an  important  mechanism  in  these  environments. 
However,  none  of  the  C„H2,,+i+  ( n  =  1—5)  alkyl  ions  exhibited 
any  appreciable  reactivity  with  O2.  All  of  these  cations  are 
closed-shell  singlet  molecules,  and  O2  is  triplet  in  its  ground 
state.  The  production  of  the  OH  radical  is  spin-allowed  and  is 
exothermic  by  over  90  kJ/mol  for  the  C„H2b+i+  (n  =  1  -4)  alkyl 
ions.  For  example, 

j-C3H7+('A)  +  02(32)  —  CH2C(OH)CH3+(2A")  +  OH(2ri) 

(7) 

is  spin-allowed  and  exothermic  by  99  kJ.  The  reaction  products 
listed  in  reaction  7  require  significant  rearrangement,  possibly 
though  a  strongly  attractive  intermediate  complex.  Several  other 
exothermic  reaction  channels,  requiring  less  rearrangement,  are 
available  but  are  spin-forbidden,  namely, 

j-C3H7+  +  02  —  CH30+  +  C2H40,  AH  =  -267  kJ/mol 
'  ‘  (8a) 

—  C2H30+  +  CH3OH,  AH  =  -249  kJ/mol 

(8b) 

—  HC0+  +  C2H60,  AH  =-209  kJ/mol 

(8c) 

It  is  not  necessary  that  a  reaction  proceed  in  a  manner  that 
allows  for  spin  conservation.  There  are  several  ion-molecule 
reactions  that  have  been  observed  to  be  fast  despite  violating 
spin  conservation.  Examples  include 

oYs)  +  co/2) — o2Yn)  +  co(‘z)  (9) 
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which  proceeds  at  the  collision  rate,49-50  the  reaction  of  CH3+ 
with  oxygen  atoms,1 ’4,? 

CH3YA)  +  0(3p)  —  HCOYA)  +  H/2)  (10) 

with  a  measured  rate  constant  of  4.4  x  10~10  cm3  s-1,  a  number 
of  positive-ion  charge-transfer  reactions  where  the  rate  constant 
exceeds  the  spin-conserving  weighted  collision  rate  constant,51 
as  well  as  a  number  of  hydrocarbon  ion  reactions  with  N(4S) 
atoms.52 

Ab  initio  calculations  regarding  C2H5+  and  j-C3H7+  reacting 
with  O2  at  the  B3LYP/6-31G(d)  level  of  theory  set  suggest  that 
these  alkyl  ions  yield  relatively  weak  complexes  with  02  with 
elongated  C-0  bond  distances,  namely, 

c2h5Ya)  +  0/2)  -  [OOC2H5]+*(3A"), 

AH  =  -24  kJ/mol  (11) 
and 

j-C3H7YA)  +  0/2)  —  [OOCH(CH3)CH3]+*(3A"), 

AH=  -13  kJ/mol  (12) 

The  minimum  calculated  C-0  bond  distance  for  these  com¬ 
plexes  is  ca.  2.6  A  with  further  C— O  approach  being  repulsive, 
independent  of  the  orientation  of  02  relative  to  the  cation.  These 
calculations  suggest  that  no  C-0  bond  can  be  formed  on  the 
triplet  surface  in  these  complexes,  and  that  the  singlet  state 
required  for  bond  formation  is  about  109  kJ/mol  higher  in 
energy.  In  addition,  hydrogen  atom  and  proton  abstraction, 
which  are  spin-allowed,  are  not  at  all  favorable  processes  due 
to  large  endothermicities.  For  example,  the  reactions  involving 
the  ethyl  cation  are 

c2h5Ya)  +  0/2)  -  OOH(2A")  +  c2h4  Yb3u), 

A H=  166  kJ/mol  (13a) 
and 

C2HsYA)  +  0/2)  -  OOHYa")  +  C2H4('Ag). 

AH  ~  248  kJ/mol  (13b) 

Therefore,  no  appreciable  O2  oxidation  chemistry  can  be 
foreseen  at  300  K  for  these  species  because  our  quantum 
chemical  calculations  predict  that  the  approach  of  the  triplet 
oxygen  molecule  toward  the  hydrocarbon  cation  is  strongly 
repulsive  and  that  the  singlet/triplet  crossing  is  not  likely  to  be 
an  important  factor  at  thermal  energy.  This  theoretical  picture 
is  consistent  with  the  present  experimental  observations,  which 
places  an  upper  bound  to  the  reactivity  of  O2  with  these 
hydrocarbon  ions  at  <5  x  10-13  cm3  s”1.  Scott  et  al.7  observed 
that  C2H5+  also  did  not  react  with  triplet  atomic  oxygen,  even 
though  the  product  channel  HCO4-  +  CH4  is  exothermic  by  399 
kJ/mol .  These  authors  suggested  that  the  barriers  on  the  potential 
energy  surface  are  too  high  for  the  extensive  rearrangement 
required.  Ozone,  on  the  other  hand,  is  able  to  produce  these 
products  in  its  reaction  with  C2H5+,  where  this  channel  is 
exothermic  by  294  kJ/mol  (Table  1). 

Ozone  is  a  singlet  in  its  ground  state,  and  ab  initio  calculations 
at  the  B3LYP/6-31G(d)  level  of  theory  indicate  that  reactant 
complexes  can  be  formed,  decreasing  in  stability  with  the  size 
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-54.5  kJ/mot 


CO,(CH3)/  !a 
-10.0  U/mo] 


-00.2  kJ/mol 


Figure  5.  Initial  reactant  complexes  for  alkyl  cations  with  ozone.  Geometries  were  optimized  at  the  B3LYP/6-3  lG(d)  level  of  theory,  bond  distances 
are  in  A,  and  bond  angles  are  in  degrees.  Interaction  energies  are  in  kJ  mol"1  with  respect  to  isolated  alkyl  cations  and  ground-state  ozone. 


R,CO+  +  RjRj 
RiCCr  +  RjR, 
RaCO^RjRj 


Figure  6.  Proposed  reaction  pathway  for  O-atom-transfer  reaction  from  ozone  to  alkyl  ions.  The  groups  labeled  Ri  to  R3  (R  =  H,  CH3,  C2H3) 
connote  the  species  surrounding  the  carbon  atom  undergoing  nucleophilic  attack  by  ozone.  After  nucleophilic  addition  of  03,  rearrangement  to  a 
transition-state  structure  enabling  02  loss  follows,  yielding  an  ionized  alkyl  oxide.  From  this  structure,  rearrangement  can  occur,  leading  to  the 
product  combinations  listed. 


of  the  alkyl  chains  attached  to  the  cationic  carbon  atom.  The 
complexes  are 

CH3+('Alg)  +  03(‘A,)  -  [CH303]+*(,A'), 

A H=  —248  kJ/mol  (14a) 
C2H5+(1A,)  +  OjOAj)  -  [CH3CH203]+(iA,)> 

AH  =  —144  kJ/mol 

(14b) 

j-C3H7+(‘A)  +  03('A,)  —  [CH3CH03CH3]+*('A), 

AH  =  -60  kJ/mol  4 
(14c) 

s-  C4H9+('A)  +  03(‘A,)  —  [CH3CH03C2H5]+*('A), 

AH  =  —55  kJ/mol 

(14d) 

r-C4H9+('A)  +  03('A,)  -  [C03(CH3)3]+*('A), 

AH  =-10  kJ/mol 

(14e) 

All  of  these  systems  form  a  complex  with  a  C— O  distance  of 
about  1.5  A  with  no  activation  energy  required  (see  Figure  5). 
From  these  complexes,  rearrangement  can  occur,  leading  to 
several  different  spin-conserving  products.  In  addition,  the  large 
exothermicity  of  the  available  reaction  channels  permits  breaking 
of  C—WC—C  bonds,  favoring  the  production  of  smaller 
fragments,  consistent  with  the  experimental  observations.  These 
complexes  are  characterized  by  weak  O3— /1-hydrogen  interac¬ 
tions;  a-hydrogens  are  not  in  the  position  to  interact  with  the 
terminal  oxygen  atom  in  the  initial  complexes.  The  presence 


of  a-hydrogens  is  therefore  not  a  factor  in  the  relative  stabilities 
of  these  complexes.  Drastic  energetic  destabilization  of  these 
initial  complexes  occurs  when  additional  alkyl  groups  are  added 
to  the  cationic  carbon  atom,  leaving  the  terf-butyl  cation/ozone 
complex  with  only  a  small  binding  energy  of  10  kJ/mol. 

A  major  product  for  all  of  the  primary  and  secondary  alkyl 
cations  is  HCO+.  The  reaction  pathway  for  this  product  channel 
could  be  traced  computationally  for  the  methyl  cation.  In  the 
case  of  CH3+,  a  straightforward  mechanism  evolved  involving 
nucleophilic  addition  of  ozone  to  CH3+,  early  0—02  cleavage 
in  the  transition  state  with  a  low  barrier  of  only  50  kJ/mol,  and 
subsequent  hydrogen  elimination  to  yield  the  highly  exothermic 
products  HCO+  +  H2  +  O2.  The  calculations  indicate  that  H202 
is  not  formed  in  this  channel  and  that  the  ethyl,  j^c-propyl,  and 
sec-butyl  cations  proceed  via  a  similar  mechanism,  as  illustrated 
in  Figure  6.  For  C2H5+,  HCO+  and  C2H30+  are  produced  from 
the  elimination  of  CH4  and  H2,  respectively,  from  the  [OCH2- 
CH3]+*  intermediate.  Similarly,  HCO+  and  C2H30+  are  pro¬ 
duced  from  the  elimination  of  C2H6  and  CH4,  respectively,  from 
the  [OCH(CH3)CH3] + *  intermediate.  On  the  basis  of  the  same 
mechanism,  the  [OC(CH3)2CH3]+*  intermediate  produced  in  the 
reaction  of  s-C4H9+  is  expected  to  produce  HCO+,  C2H30+, 
and  C3H50+  with  the  elimination  of  C3Hg,  C2H6  and  CH4, 
respectively.  This  mechanism  is  in  qualitative  agreement  with 
the  branching  fractions  observed  in  these  reactions,  as  shown 
in  Table  1,  although  only  a  barely  detectable  amount  of  C3H50+ 
reaction  product  is  observed  in  the  scc-butyl  cation  reaction  with 
ozone,  perhaps  due  to  steric  constraints.  This  mechanism  is 
similar  to  that  proposed  by  Mendes  et  al.11  on  the  basis  of  their 
ab  initio  calculations  of  O-atom-transfer  reactions  of  ozone  with 
radical  cations  of  pryridines,  alkyl  halides,  and  halogen  cat¬ 
ions. 


1988  J.  Phys .  Chem .  A,  Vol.  108,  No.  11,  2004 


Williams  et  al. 


A  possible  explanation  for  the  apparent  lack  of  reactivity  of 
t-C^Ug*  and  r-CsHu*  is  that  the  initial  complex  with  ozone  is 
not  formed.  The  calculations  find  that  the  stability  of  the 
complex  decreases  from  248  to  10  kJ/mol  on  going  from  CH3+ 
to  f-C4H9+.  Therefore,  it  seems  likely  that  the  lack  of  reactivity 
of  the  alkyl  ions  involving  a  tertiary  carbocation  may  be 
attributed  to  a  greatly  reduced  O3  binding  energy. 

V.  Summary 

Numerous  alkyl  cations  and  HCO+  reacting  with  O2  and  O3 
were  studied  at  300  K  in  a  selected-ion  flow  tube  coupled  to  a 
novel  ozone  source.  The  HCO+  cation  exhibits  no  reactivity 
with  O2,  but  the  reaction  of  HCO+  with  O3  proceeds  at  the 
collision  rate  via  proton  transfer  to  produce  H03+.  None  of  the 
alkyl  carbocations  studied  showed  any  measurable  bimolecular 
reactivity,  k  <  5  x  10”13  cm3  s-1,  with  O2,  despite  the 
availability  of  reaction  channels  with  exothermicities  of  several 
hundred  kilojoules  per  mole.  The  ab  initio  calculations  indicate 
that  the  O2  reaction  systems  form  weak  complexes  with  large 
C— O  bond  distances  on  the  lowest  energy  triplet  potential 
energy  surface  which  are  repulsive  at  smaller  distances.  Access 
to  the  singlet  surface  is  required  for  bond  formation;  however, 
this  surface  is  not  accessible  at  thermal  energies. 

For  reactions  with  O3,  the  total  reaction  rate  constants  for 
the  alkyl  carbocations  decrease  dramatically  as  the  order  of  the 
reactant  carbocation  increases,  as  shown  in  Figure  2  and  listed 
in  Table  1 .  This  finding  is  consistent  with  ab  initio  calculations 
which  show  that  the  initial  alkyl  cation/ozone  complexes  are 
largely  destabilized  when  alkyl  groups  are  added  to  the  cationic 
center.  The  primary  carbocations  CH3+  and  C2H5+  react  at  100% 
and  46%  of  the  collision  rate  given  by  the  Su— Chesnavich 
equation,  respectively.  The  secondary  .y-CaH?4*  and  $-C4H9+ 
carbocations  reacted  with  O3  at  19%  and  7%  of  the  collision 
rate,  respectively.  The  tertiary  carbocations  f-C4H9+  and  f-C* Hn+ 
were  found  to  be  unreactive  with  O3,  k  <  5  x  10“ 12  cm3  s_1, 
which  is  the  detection  limit  of  our  apparatus  using  this  ozone 
source.  A  straightforward  mechanism  was  discussed  involving 
initial  complex  formation  involving  the  primary  or  secondary 
alkyl  cation  with  O3,  followed  by  early  0-02  cleavage  in  the 
transition  state,  and  finally  hydrogen  and  alkane  elimination  to 
yield  the  highly  exothermic  products  listed  in  Table  1.  An 
explanation  for  the  lack  of  reactivity  of  f-CUHk4'  and  r-C5Hn+ 
is  that  the  initial  complex  with  ozone  is  not  formed  because  of 
a  much  reduced  O3  binding  energy.  Details  of  theoretical  studies 
of  the  potential  energy  surfaces  of  the  present  systems  will  be 
published  separately. 

The  only  other  positive  organic  ion  reactions  with  ozone 
previously  studied  used  a  pentaquadrupole  mass  spectrometer 
at  nonthermal  (ca.  1  eV)  collision  energies,  and  thermal  rate 
constants  were  not  reported  as  done  here.11  In  addition,  alkyl 
ions  were  not  part  of  that  study.  Therefore,  the  values  reported 
in  Table  1  represent  the  most  self-consistent  and  complete  set 
of  reaction  rate  and  branching  fraction  data  reported  to  date 
and,  on  the  basis  of  the  above  considerations,  are  recommended 
for  use  in  modeling  applications  requiring  thermal  rate  data  near 
300  K. 
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